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We describe the transition from a ferromagnetic phase, to a disordered paramagnetic phase, which 
occurs in one-dimensional Kondo lattice models with partial conduction band filling. The transition 
is the quantum order-disorder transition of the transverse-field Ising chain, and reflects double- 
exchange ordered regions of localized spins being gradually destroyed as the coupling to the conduc- 
tion electrons is reduced. For incommensurate conduction band filling, the low-energy properties of 
the localized spins near the transition are dominated by anomalous ordered (disordered) regions of 
localized spins which survive into the paramagnetic (ferromagnetic) phase. Many interesting prop- 
erties follow, including a diverging susceptibility for a finite range of couplings into the paramagnetic 
phase. Our critical line equation, together with numerically determined transition points, are used 
to determine the effective range of the double-exchange interaction. Models considered are the spin 
1/2 Kondo lattices with ferromagnetic and antiferromagnetic couplings, and the Kondo lattice with 
repulsive interactions between the conduction electrons. 



I. INTRODUCTION 

The Kondo lattice model (KLM) describes the interac- 
tion between a conduction band and a lattice of localized 
spins. The hamiltonian for the KLM is 

H = -t J2 E c U CT + J E s ^- s ^ (!) 

<i,j> a j 

where t > is the conduction electron (or simply elec- 
tron) hopping, and < i,j > denotes nearest neighbors. 
Sj are spin 1/2 operators for the localized spins, and 

^cj = a ' c \a CT < J - l7 ' C 3tj' are pseudospin operators for 

the conduction electrons, er arc Pauli spin matrices, and 
Cj a ,Cj a the electron site operators. We choose units so 
that the hopping t = 1, and measure the Kondo cou- 
pling J in units of t. The KLM is an effective model 
for heavy-fermion systems when the coupling J is small 
and positive and models colossal magnetoresistance 
(CMR) materials when J is large and negative |Q. 

In the following, we give a comprehensive description of 
the transition from a ferromagnetic (FM) ordering of the 
localized spins at stronger couplings, to a quantum dis- 
ordered paramagnetic (PM) phase at weaker couplings. 
The transition has been identified in the onc-dimcnsional 
(ID) KLM with partial conduction band filling both an- 
alytically |||, and in numerical simulations by a variety 
of methods We focus mainly on the J > ID 

KLM relevant to heavy-fermion systems, in which the lo- 
calized spins model /-electrons in lanthanide or actinide 
compounds. The bulk of the numerical simulations Q-j?]] 
are devoted to this model. Partial conduction band fill- 
ing n — N c /N < 1 is assumed throughout, where TV C 
is the number of electrons, and TV the number of local- 
ized spins. Our approach generalizes to ID KLMs with 
repulsive interactions between the electrons, and also to 



the KLM with a FM coupling J < 0. These models are 
considered in Sec. IV. 

A summary of our method is as follows: Abelian 
bosonization is used to describe the conduction band. 
Using a unitary transformation, the bosonized KLM is 
written in terms of a basis of states in which the localized 
spin and electron spin degrees of freedom are coupled. In 
the new basis, the competing interactions leading to the 
FM-PM transition are clearly exhibited. The compet- 
ing effects are double-exchange FM ordering at stronger 
coupling, and spin-flip disorder processes at weaker cou- 
pling. Since the FM-PM transition is signalled by the 
ordering of the localized spins, we take expectation val- 
ues for the electron Bose fields, and obtain an effective 
hamiltonian for the localized spins. The effective hamil- 
tonian maps to a transverse-field Ising chain close to the 
FM-PM phase boundary, and we determine the critical 
line for the resulting quantum order-disorder transition, 
as well as many properties of the localized spins near crit- 
icality. At weak-coupling deep in the PM phase, the effec- 
tive hamiltonian determines RKKY-like behavior, with 
dominant correlations in the localized spins at 2/cp of 
the conduction band. 

The new element here is an emphasis on the double- 
exchange interaction, which tends to align the localized 
spins at stronger couplings. Double-exchange is often ig- 
nored in discussions of the J > KLM. Following its 
introduction 20 years ago by Doniach ||, the J > 
KLM has usually been discussed in terms of a compe- 
tition between Kondo singlet formation, and the RKKY 
interaction. For the KLM with a half-filled conduction 
band n = 1, this characterization appears sufficient. For 
the KLM with a partially-filled conduction band n < 1, 
it has been realised that neither Kondo singlet forma- 
tion, nor the RKKY interaction, are sufficient to describe 
rigorously established properties of the model [pi. A 



1 



succession of analytic |lll|l2l[| and numerical re- 
sults have established an extensive region of FM order- 
ing at stronger coupling. This cannot be explained in 
terms of RKKY, which operates at weak-coupling, nor 
in terms of Kondo singlets, since they are non-magnetic. 
The missing element is double-exchange ordering due to 
an excess of localized spins over conduction electrons. 
Double-exchange requires only that N > N c (i.e. n < 1), 
and that the conduction electron hopping t ^ 0. It op- 
erates in any dimension, for any sign of the coupling J, 
and for any magnitude Sj of the localized spins Q . The 
double-exchange interaction is specific to the Kondo lat- 
tice, and is absent in single- or dilute-impurity systems in 
which the situation is reversed, and the electrons greatly 
outnumber the localized spins. 

Double-exchange is conceptually a very simple inter- 
action: For n < 1 each electron has on average more 
than one localized spin to screen, and consequently hops 
between several adjacent spins gaining screening energy 
at each site, together with a gain in kinetic energy. 
Since hopping is energetically most favorable for elec- 
trons which preserve their spin as they hop (called coher- 
ent hopping) , this tends to align the underlying localized 
spins jn|. We show in Sec. II that coherent conduction 
electron hopping over a characteristic length a may be 
incorporated into a bosonization description which keeps 
the electrons finitely delocalized. At lengths beyond a, 
the electrons are described by collective density fluctu- 
ations, as is usual in ID Fermi systems. The electrons 
remain finitely delocalized over shorter lengths, and de- 
scribe coherent hopping over several adjacent sites. This 
tends to align the underlying localized spins at stronger 
coupling, a measures the effective range of the double- 
exchange interaction, and is in principle a function of 
both filling n and coupling J. 

Our approach generates a ground-state phase diagram 
for the partially-filled ID KLM in agreement with avail- 
able exact and numerical results for the J > KLM 
Q^IJ] §, for the J > KLM with repulsive interac- 
tions between the conduction electrons JlJ,|l^|, and for 
the KLM with a FM coupling J < f|. Interesting new 
properties we determine include the following: 

(i) For incommensurate conduction band filling, anoma- 
lous regions of double-exchange ordered localized spins 
survive close to the transition in the PM phase. Simi- 
larly, anomalous regions of disorder survive close to the 
transition in the FM phase. These regions exist due to 
the inability of the incommensurate conduction band to 
either totally order, or totally disorder the localized spins 
as the transition is crossed. Although the anomalous re- 
gions are very dilute, they dominate the low-energy prop- 
erties of the localized spins. Many interesting results fol- 
low, including a diverging susceptibility in the PM phase 
for a finite range of couplings close to the transition. 

(ii) The effective range of double-exchange ordering, mea- 
sured by the length a for coherent electron hopping, is 



determined using our critical line equation together with 
numerically determined transition points. We find that 
a scales as 1/ y/\J\ for any sign of the coupling J on the 
transition line. For J < 0, a — > as half-filling is ap- 
proached; double-exchange is ineffective as n — * 1, and 
the transition line diverges. For J > 0, a is reduced to 
less than a lattice spacing as n — ► 1, but remains non- 
zero; the transition line remains finite close to half-filling 
for the J > KLM. 

An outline of the paper is as follows: In Sec. II we de- 
rive the effective hamiltonian for the localized spins in the 
partially-filled J > ID KLM, and provide a justification 
of results briefly reported elsewhere H . In Sec. Ill the ef- 
fective hamiltonian is analyzed to determine the ground- 
state magnetic phase diagram, together with low-energy 
properties of the localized spins near the resulting FM- 
PM transition. Using our results, together with available 
numerical data, we determine the effective range of the 
double-exchange interaction on the transition line. Our 
method may be generalized to describe the FM-PM tran- 
sition in other KLMs, and in Sec. IV we consider the ID 
KLM with interactions between the electrons, and show 
that the FM-PM transition is pushed to lower values of 
J for repulsive interactions. We consider also the KLM 
with a FM J < coupling, and follow a similar analy- 
sis to that for the J > case. The ground-state phase 
diagram is determined, together with the effective range 
of the double-exchange interaction on the resulting FM- 
PM transition line. The paper concludes in Sec. V with 
a summary of results, and a discussion of the general fea- 
tures of the FM-PM transition in partially-filled KLMs. 

II. EFFECTIVE HAMILTONIAN FOR THE 
LOCALIZED SPINS 

A large class of ID many-electron systems may be 
described using bosonization techniques Jl6|: The elec- 
tron fields may be represented in terms of collective den- 
sity operators which satisfy bosonic commutation rela- 
tions. Bose representations provides a non-perturbative 
description which, in general, is far easier to manipulate 
than a formulation in terms of fermionic operators. In 
the ID KLM, the conduction band may be bosonized, 
but not the localized spins. This is because the spins 
are strictly localized, and their Fermi velocity vanishes. 
Moreover, since there is no direct interaction between 
the localized spins in the KLM, it not possible to use 
bosonization via a direct Jordan- Wigner transformation. 
We therefore bosonize only the conduction band. 

Bose representations are conveniently written in terms 
of Bose fields, defined as follows: 

Mj) = -*E Tf [v+{k) + M*)]M*)e <fcio > 
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(2) 



In Eqs. (|2j) v — p,(j labels charge and spin, with charge 
and spin density operators 



r (k) = 2J Pra(k), U r {k) 



ap ra (k). 



(3) 



(j = ±1 as the spin is up or down, respectively. The 
density operators p rC r(k) are the basic bosonic objects, 
and are defined by 



Praik) 



0<rfc'<7r/a 



The density operators describe collective coherent 
particle-hole excitations about the right (r = +) and 
left (r = — ) Fermi points at +kp and — hp, respectively. 
(hp = im/2a with a the lattice spacing. The system 
length is L = Na.) The density operators p rC r(k) are 
bosonic for wave vectors |fc| up to a -1 . For these wave 
vectors we have (T^l 



[Pra{k),Pr'a'{k')] = 5 rr >S k k > 5 a y . (5) 

In bosonization, a measures the minimum wavelength for 
the densities p rr j(k) which satisfy the bosonic commuta- 
tion relations Eq. (^). Straightforward calculation, as in 
Ref. [17], shows that a must satisfy a > O(kp)^ 1 . This 
is clear also on physical grounds; fermionic density op- 
erators are not collective, and hence cannot be bosonic, 
at wavelengths of the order of the average interparticle 
spacing. In the Bose fields of Eqs. (En, A a (k) is a cut-off 
function on bosonic density operators. The cut-off func- 
tion is an even function of k, and satisfies A a (fc) « 1 when 
\k\ < a^ 1 , and A a (k) when \k\ > a^ 1 . A a (k) ensures 
that only bosonic density operators enter the Bose fields. 
The commutation relations between the Bose fields are 
then c-numbers. The physical significance of the Bose 
fields is as potentials: d x <p v (j) is proportional to the v- 
density at site j, and d x 9 u (j) = H v (j) is proportional 
to the average z/-current at j. (d x ip v (j), ip = (f),9, is 
shorthand for d x il) u {x/a) evaluated at x = ja.) 

Bose representations for Fermi operators are derived 
by requiring that they correctly reproduce the commu- 
tation relations of the Fermi operators with the density 
operators p ra -(k), and by requiring that they correctly 
reproduce the non-interacting (J = 0) expectation val- 
ues. Since the states generated by the density operators 
Pra(k) span the ID state space |18 , this prescription en- 
sures that the Bose representations will reproduce the 
same matrix elements as the original Fermi operators. In 
this way, the hopping term in Eq. (jl]) is given by 
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^2{nl(j) + [d x Mj)} 2 } 



(6) 



I'-J 



with a linearized dispersion, and Fermi velocity vp = 
2asin(7m/2) in units of t. To bosonize the KLM, we 
also require representations for the on-site Fermi bilin- 
J2 r .r' c l-jcr c r>io' ■ The off-diagonal bilinears 



ears Cj a c j<7 
J. «. , • 

-"rja r'jer 



,, in which r r' and/or a ^ a', may be con- 
structed from the Bose representation for the single Fermi 
site operators 



Crja- = y Aa/2a expi[rkpja + ^ r a(j)], 

*«r0') = {°pU) + r Mj) + <r[9*U) + r^0')]}/2, (7) 

where A is a dimensionless constant depending on the 
cut-off function A a (k). The representation Eq. (^) for 
(4) c r j a reproduces the correct commutation relations only 
with long- wavelength density operators p ra (k). Thus the 
representation Eq. (^) may not correctly reproduce the 
short-range properties of the original conduction electron 
site operators Cj a . The Bose representation for the diag- 
onal on-site bilinears, i.e. the density operators in real 
space, may be obtained directly from their Fourier ex- 
pansion: 



Vc^ c ■ 



= ^d x [(j) p (j) + act 



r(j)] 



(8) 



to an additive constant depending on n. As for Hq, and in 
constrast to Crja, this representation is exact. Substitut- 
ing these representations into Eq. ([[]) gives the bosonized 
KLM hamiltonian 
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A ^ { cos [^(i)] + cos[2k F ja + <f> p (j)]} 



Ja 



x (e- i6 ^Sj +h.c.) 



- A— V sin[^(j)] sm[2kpja + cf> p (j)]S] . (9) 
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The bosonized hamiltonian generates the same behav- 
ior as the original hamiltonian provided that the conduc- 
tion electrons are not strongly localized. In particular, 
the bosonized hamiltonian does not directly describe res- 
onant on-site Kondo scattering: At strong-coupling the 
electrons localize, with each forming a Kondo singlet with 
the localized spin at the same site 19,12]. The localized 
singlet formation is governed by the short-range proper- 
ties of the spin-flip bilinears Cj a Cj_ a . However, the Bose 
representation for these terms, derived from Eq. (j^), is 
reliable only at long- wavelengths. It describes the prop- 
erties of a spin-flipped electron only at large distances 
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from the scattering site. This provides a good descrip- 
tion at weaker couplings, as usual in the bosonization of 
ID Fermi systems, but may be insufficient when the cou- 
pling is strong enough that the electron becomes trapped 
on-site by the localized spin. 

A second point to note about the bosonized hamilto- 
nian concerns spin-rotation symmetry. SU (2) symmetry 
is manifest in the original hamiltonian, Eq. ([!]), for both 
the electrons and the localized spins, but is obscured in 
the bosonized version. This is due to the use of abelian 
bosonization, which treats the electron spin z direction 
on a special footing, and breaks the SU(2) electron spin- 
rotation symmetry down to £7(1) . To see the effect of 
this, note that the original hamiltonian preserves both 
the total spin S to t as well as its z component S* ot , and 
at stronger couplings in the FM phase may be decoupled 
into subspaces with different values of 5* ot . (See, for 
example, Ref. [12].) Abelian bosonization effectively sin- 
gles out the subspace with maximal Sfot in the FM phase 
(cf. Eq. ([l^) below), and may be physically motivated by 
crystal-field effects. 



Using [Sf, S?,] = iSjj'Sj etc., we get 



so that S rotates the localized spins in the xy-plane de- 
pending on the local electron spin current field. S does 
not alter the Sj configuration. The only Bose field trans- 
formed under S is the electron spin density field, which 
takes a form that depends on the localized spin configu- 
ration all along the chain: 

Mj) = Mj) + K(j), 

where K(j) is a long-range object which essentially 
counts all the Sp to the right of site j, and subtracts 
from that the S% to the left of j: 

K(j) = i^Yy<t>„{jie a (j>)]s* r . (ii) 



A. Unitary transformation 

A simple method for determining the ordering induced 
on the localized spins by the electrons is to choose a basis 
of states in which competing effects become more trans- 
parent. This is achieved by applying a unitary transfor- 
mation which changes to a basis of states in which the 
conduction electron spin degrees of freedom are coupled 
directly to the localized spins. We choose the transfor- 
mation 

exp(S), S = i^£>(j)SJ. 

3 

A variant of this transformation was first used by Emery 
and Kivelson for the single- impurity Kondo problem [po[ , 
and was later gener alized to the ID KLM [|IJ. The us- 
age here is different. Ref. [21] aimed to describe the 
conduction electrons, and the transformation was used 
to remove the spin current field 6 a (j) from the hamilto- 
nian. Here we aim to describe interactions between the 
localized spins, in which the electrons act as the media- 
tors. The form of the transformation is then chosen so 
as to make explicit a FM ordering of the localized spins. 
This effect was entirely missed in the previous work [ pl| . 
The factor Ja/2nvF is chosen so that terms of the form 
[dx<l>i?(j)]Sj exactly cancel in the transformed hamilto- 
nian. This permits the ground-state Sj configuration to 
be chosen independent of the on-site electron spin density 
in the transformed basis. 

Transformed operators O — e~ s Oe s may be calculated 
using the standard commutator expansion 



Here the Bose field commutator [4> a (j), 0o{j')\ sign(j — 
j 1 ) in for (j — j')a ^> a, and this provides the main con- 
tribution to K{j). This term is discussed further in Sec. 
III. Related to the transformation of the spin density field 
is the transformation of the actual electron spin density 
at j. This depends on the local Sj, configuration. The 

transformed spin density d x (j) a (j) is given by 



Ja 

TTVf 



dk cos[k(j-j')a]A 2 a (k)}Sp, 



The integral here is the Bose field commutator 
i[</> a (J),H <r (J')]/2, and is discussed further below. 

After some manipulation, the above results give the 
transformed KLM hamiltonian 
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J z a 



Ait 2 vf 

3,3' 

Ja ■ 



dk cos[k(j - f)a]Al(k) }S*Sp 



o 



= + [0,S]/l! + [[0,S],S]/2! 



(10) 



+ A ^ i cos [^(i) + Mi)] + cos[2k F ja + p (j)]} 

3 

x (e-^+^/^^St + hx.) 
- A— V sinTO) + Mj)\ sxn[2k F ja + Mj)\S- 

(12) 

provided that the cut-off function is not too 'soft': 
A™(fc) ps A a (k),m — 2,3,4. (Discrepancies near |fc| w 
a -1 introduce negligible corrections.) Note that the uni- 
tary transformation has been carried out exactly and not 
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perturbatively, i.e. there has been no artificial truncation 
of the commutator series of Eq. (|l^). (The c-number 
Bose field commutators are essential for this.) It follows 
that the transformed hamiltonian of Eq. ([jj]) is identical 
to the bosonized hamiltonian, and is an exact rewriting 
of Eq. (||) in terms of a new basis of states in which the 
conduction band and localized spins are interwoven. 



B. Double-exchange ordering 

The important new term in the transformed hamilto- 
nian Eq. (|l^) is the second: 



2^2 



J 2 a 
4tt 2 vf 



E 

3,3' 



dk cos[k( 3 - j')a]A 2 a {k) 



QZ QZ 

°3 °3' 



(13) 



It represents a non-perturbative effective interaction be- 
tween the localized spins, and is the only one of this 
type to be derived for the KLM; other effective interac- 
tions, namely the RKKY interaction at weak-coupling, 
and the strong-coupling effective interaction of Sigrist et 
al. Q, are both perturbative. We consider FM in the 
KLM in some detail in this subsection. First we analyze 
the properties of the interaction described by Eq. (13), 



and describe how it arises from the bosonization of the 
conduction band. Second we present previously known 
properties of the double-exchange interaction. The inter- 
action of Eq. ([l3|) shares these properties, and we identify 
it as the double-exchange interaction in the KLM. Since 
double-exchange is usually not considered in discussions 
of the KLM with J > 0, we conclude the subsection with 
a simple intuitive picture of double-exchange ordering in 
the KLM at low conduction band filling. 

Eq. ( |l3| ) possesses the following properties: 
(i) the term originates, via bosonization and then the uni- 
tary transformation, from the terms Hq and the forward 



scattering part of (J/2) J2j( n j 1 i — n jl)^j ( n . 



in the KLM hamiltonian Eq. (|lj). (Note that the Bose 
representations for the electrons in these terms are ex- 
act.) 

(ii) Eq. ( p"3] ) is independent of the sign of J, and takes the 
same form for any magnitude S 2 of the localized spins. 

(iii) Since Eq. ( |l3| ) is of order J 2 , whereas the remaining 
terms in the transformed hamiltonian Eq. (|l|) are of or- 
der J, the interaction Eq. ( p"3] ) dominates the ordering of 
the localized spins as J increases. 

(iv) Eq. (^||) is FM for all (differentiable) choices of the 
cut-off function A Q (fc). 

To give examples of the form of the FM interaction Eq. 
( |l3| ) in real space, consider Gaussian and exponential cut- 
off functions defined by A a (k) = exp(— a 2 k 2 /2) for a 
Gaussian cut-off, and by A a (k) = exp(—a\k\/2) for an 
exponential cut-off. For these cut-off functions, the inte- 
gral in Eq. (O) reduces to 



dk cos(fc ja) A 2 (fc) = (y/n/2a) exp ~(ja/2a) 2 (14) 

for the Gaussian, and 

> 

dk cos(kja)A 2 a {k) = a/ (a 2 + (ja) 2 ) (15) 



for the exponential. The integrals are positive and non- 
negligible for ja < a. The form of the FM interaction 
for Gaussian and exponential cut-off functions is shown 
in Fig. 1. It is clear that the length a characterizes the 



effective range of the FM interaction of Eq. (13) 



Gaussian cut-off 

v\ Exponential cut-off 




FIG. 1. The range in real space of the FM interac- 
tion Eq. (113) for exponential exp— (a|fc|/2) and Gaussian 
exp — (oP'krf'l) cut-off functions A a (k). J c g is the interaction 



) strength in units of aj a /4-7T vf- 



The interaction Eq. ( |l3|) originates from the bosoniza- 
tion of the conduction band as follows: At wavelengths 
beyond a, the electrons are involved in collective den- 
sity fluctuations Eq. (|J) . These fluctuations involve large 
numbers of electrons, and satisfy bosonic commutation 
relations Eq. (||). This is the standard behavior of ID 
many-electron systems for weak interactions. At wave- 
lengths below a, the density fluctuations are not collec- 
tive, and do not satisfy bosonic commutation relations. 
Since bosonization describes fluctuations only over sepa- 
rations beyond a, the bosonization description is equiv- 
alent to keeping the electrons finitely delocalized over 
a, with the electrons preserving their spin over this 
range. Eq. ( |l3| ) is the ordering consequently induced 
on the localized spins by the finitely delocalized elec- 
trons, and arises formally from the Bose field commu- 
tator [4> a (j), H a (j')]. This commutator takes canonical 
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(5-function form in field theory []16|, but is smeared over 
a range a > O(kp)^ 1 for the conduction band, which 
has a finite density n of electrons. The smearing reflects 
the inability of the Bose fields to distinguish separations 
below a. 

We turn now to briefly summarize previously known 
properties of the double-exchange interaction. Double- 
exchange was first proposed many years ago by Zener 
PH to explain FM ordering in mixed-valency mangan- 
ites, and is of much current interest in relation to CMR 
materials The essential characteristic required for a 
system to exhibit double-exchange ordering is that the 
number of electrons be less than the number of local- 
ized spins. In this case, consider first infinitely strong 
coupling J — oo. Each electron forms a perfectly local- 
ized on-site spin singlet (or triplet for J < 0) with the 
localized spin at the same site. The remaining N — N c 
unpaired localized spins are free. When the conduction 
electron hopping is turned on, the electrons gain energy 
by hopping to unoccupied sites, since they gain energy 
both by screening the unpaired localized spin, together 
with a gain in kinetic energy. Since electrons tend to 
preserve their spin as they hop, called coherent hopping 
Pj , this tends to align the underlying localized spins [jl3| . 
This is the double-exchange mechanism, and is generated 
by the conduction electron kinetic energy and the diago- 
nal part of the on-site interaction between electrons and 
localized spins. Double-exchange is always FM, and dom- 
inates at stronger couplings p3] . It is the physical basis of 
the FM rigorously established in the ID KLM by Sigrist 
et al. [ pl]p2| . (See also Ref. [10].) Since double-exchange 
ordering requires only TV > N c and a non- vanishing hop- 
ping, its existence (as opposed to other properties such 
as its effective range) does not depend on the sign of J, 
nor on the magnitude of the localized spins. This is clear 
from the early analysis of Anderson and Hasegawa [fL3t ; 
their result does not depend on the sign of J, except 
for some numerical prefactors, and is semiclassical with 
nearly trivial quantum modifications. 

Properties (i)-(iv) above for the interaction of Eq. (|T^ ) 
are identical to those of a double-exchange interaction. 
This leads us to identify Eq. ([l^) as the double-exchange 
interaction in the partially-filled ID KLM. Coherent 
conduction electron hopping, which generates double- 
exchange ordering, is described in the bosonization by 
electrons finitely delocalized over lengths a, and a mea- 
sures the effective range of the double-exchange interac- 
tion, as in Fig. 1. Note that a enters the bosonization 
description as an undetermined but finite length; from 
bosonization, we know only that a > 0(kp)^ 1 , and that 
in general a will be a function both of filling n and cou- 
pling J 0. In Sees. Ill and IV (cf. Fig. 4), we deter- 
mine a in a special case by using our results together 
with available numerical data. 

Since many discussions of the J > KLM neglect the 
double-exchange interaction, it is useful to present the 



following simple characterization valid at low conduction 
band filling. Double-exchange may be described by the 
KLM hamiltonian Eq. (|l|) with spin-flip interactions ig- 
nored: 
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The occupation of a site by an electron with the same 
spin as the localized spin costs an energy J/2. We exclude 
these states as a first approximation valid at stronger cou- 
plings. At small filling, we consider a finitely delocalized 
electron of spin a spread over sites j for which the local- 
ized spins Sj have spin —a for J > 0. From ifdc, the 
wavefunction ^> a {x) for the electron, in the continuum 
limit, satisfies the nonlinear Schrodinger equation 

d 2 x i> a {x) + (Jm/2)|V> CT (z)| 2 Vv(z) = 2mEip a (x) (16) 

with m the bare electron mass. The electron gains energy 
due to its occupation |i/v(2)| of a point with localized 
spin S* of spin —a, and this generates the non-linearity. 
Finitely delocalized solutions of Eq. (16) are solitons with 
11 

ip a (x) = B e ix sech (B^Jm/4 (x - a*)) . (17) 

B and xq are constants. Our simplified picture of the 
KLM at low conduction band filling is then of a gas of 
solitons. The solitons may be pictured as spin polarons 
0,0, and describe the dressing of each electron by a 
cloud of localized spins which align opposite to the con- 
duction electron spin for J > 0. The spatial extension 
of the polarization cloud characterizes the range of the 
indirect FM ordering induced on the localized spins by 
the electron, and is equivalent to the effective range a of 
the double-exchange interaction as described previously. 
From Eq. (|l7|), the polarization cloud decays exponen- 
tially at large distances with a characteristic length scale 
proportional to 1/vJ. This gives a low density form 
a/a oc \/\fj. At vanishingly small fillings, we may use 
the exact solution of Sigrist et al. |ll|] for the KLM with 
one conduction electron. The polarization cloud decays 
exponentially for small J in the FM phase, with a char- 
acteristic length a/a = \Jl/J. 

Since the interaction Eq. ( |l~3| ) is short-range for all fi- 
nite a (correct for all finite n, cf. Fig. 4), we approx- 
imate it in the usual way by its nearest-neighbor form 
-^Ej^+i, where 



J 2 a 2 
2tt 2 vf 



J = 



dk cos(ka)A 2 t (k) . 



(18) 



We do not expect critical properties to be affected by this 
approximation. 
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C. Effective hamiltonian 

We aim to use the transformed hamiltonian Eq. (jl^ 
to determine the ground-state properties of the local- 
ized spins in the partially-filled ID KLM. Our concen- 
tration on the localized spins is motivated by the re- 
sults of numerical simulations on the KLM. Simulations 
on large chains have been carried out on the partially- 
filled ID KLM using quantum Monte Carlo 0], and using 
the density-matrix renormalization-group (DMRG) 
The results uniformly show that the correlations between 
the localized spins are much stronger than the correla- 
tions between the conduction electrons, and the FM-PM 
transition is signalled by the crossover from FM to incom- 
mensurate (generally 2kp) correlations in the structure 
factor of the localized spins. The corresponding electron 
correlations are observed to weakly track those of the lo- 
calized spins, and the electron momentum distribution 
shows no dramatic change as the FM-PM transition line 
is crossed ](| . The freezing of the electron spin degrees of 
freedom occurs only at very strong coupling deep in the 
FM phase. 

An effective hamiltonian for the localized spins is ob- 
tained from Eq. (|l^) by taking appropriately chosen ex- 
pectation values for the conduction electron Bose fields. 
Since the Bose fields enter only in the weak-coupling 
terms of order J in Eq. (|l2|), we approximate the Bose 
fields by their non-interacting J = expectation values: 



M?)>O = <0a(?")>O = O. 



(19) 



This holds for the charge density field 4> P (j), since at 
weak-coupling the charge structure factor is free electron- 
like p|. For the spin fields, Eq. ( |l9| ) follows from real- 
space renormalization-group studies |25| ], which show 
that the spin degrees of freedom of the ID KLM flow to 
the non-interacting fixed point at weak-coupling. (This 
property is specific to the Kondo lattice p5| . For the 
single-impurity Kondo model at zero temperature, the 
spin coupling renormalizes to infinity for all J =/= 0.) 
Note that Eq. ( |l9| ) is further supported by a study of 
the ID KLM with t-J interacting conduction electrons 
p5[ . Using a combination of exact diagonalization and 
the DMRG, the same ordering was observed to be in- 
duced on the localized spins as in the pure KLM of Eq. 
([j]) , and confirms the insensitivity of the ordering to the 
details of the conduction electron behavior. The trans- 
formed hamiltonian Eq. ( |l2"|) now reduces to an effective 
hamiltonian for the localized spins: 



H c ff — ^ ' Sj Sj_^i 



Ja 



A— V{cos[if(j)] + cos[2fc F ja]}S*f 
a ^-^ ' J 



A— V sm[K(j)] sin[2k F ja]S? 



III. GROUND-STATE MAGNETIC PHASE 
DIAGRAM 

In this section the effective hamiltonian H e g of Eq. j20| ) 
is analyzed to determine the ground-state properties of 
the localized spins as a function of conduction band fill- 
ing n and Kondo coupling strength J > 0. Since the 
FM double-exchange coupling J is of order J 2 (cf. Eq. 
(|l8|)), it is immediate from Eq. ( p0| ) that H e g determines 
a FM ordering for the localized spins at stronger cou- 
plings J ^> 1 for all fillings n < 1. We find below that 
the FM ordering is gradually destroyed as the coupling 
J is lowered. The destruction of the FM order is deter- 
mined by the second term of Eq. ( |20| ) ; the effective hamil- 
tonian takes the form of a transverse-field Ising chain in 
the phase transition region, and the KLM undergoes a 
quantum FM-PM transition at a filling dependent crit- 
ical coupling J c . The critical coupling is of order unity 
at most conduction band fillings. The determination of 
the FM-PM transition, and a disussion of the properties 
of the localized spins near the phase boundary, are con- 
tained in subsection IIIA. In subsection IIIB, we consider 
the weak-coupling regime J< 1. At weak coupling the 
double-exchange ordering is ineffective, and H e g reduces 
to a system of free localized spins in fields determined 
by conduction electron scattering (the last two terms of 
Eq. (pcf)). We find that H e g determines dominant 2k f 
(RKKY-like) correlations in the localized spins at weak 
coupling. Finally, in subsection IIIC, we plot the ground- 
state phase diagram for the KLM as determined by the 
effective hamiltonian. 



A. The FM-PM phase transition 

We begin our analysis of H e s by evaluating the long- 
range object K(j) in the strong-coupling FM phase. 
K(j), given by Eq. (|ll|), originates with the unitary 
transformation S. It has some similarity to the disor- 
der term in the Jordan- Wigner transformation, but in- 
stead of counting the Sp only over sites to the left of 
j, it counts also the Sp to the right. In the thermody- 
namic limit N — > oo, and using the large j — j' form 
[4>o-{j), ^cr(j')] = sig n (i ~ j') * 7r f° r the Bose field commu- 
tator, we may rewrite Eq. (|ll|) in the form 



(21) 



/ = 1 



(20) 



where €j(l) — SJ +l — Sp_ l . has possible values 0, ±1 
for large I. For small j—f, the commutator [4> a (j), ^o-(j')] 
grows smoothly from zero at j — j', to sign(j — j') itt at 
(j — j) = 0(a/a). The exact form of the commutator 
at short-range depends on the choice of cut-off function 
Ao,(fc), but all we require here is the general form: The 
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effect of short-range corrections to the Bose held com- 
mutator is just to allow £j(l) to take values between —1 
and 1 for I < a/a. K{j) is then similarly smoothed, and 
takes values between integral multiples of Ja/2vp. 

By writing K(j) in the form of Eq. (21), it is clear that 
K(j) vanishes in the FM phase in a thermodynamically 
large system. Indeed K(j) will not be appreciable until 
the system is strongly disordered. It follows that any 
transition out of the FM phase will be governed by the 
first two terms of the effective hamiltonian Eq. (^0|) . For 
convenience we collect these terms in the hamiltonian 
ffcrit (with 'crit' for critical): 



ffcrit — 3 ^ Sj Sj_\_i 



A— Y {l + cos(2k F ja)}S* 



(22) 



ffcrit is a quantum transverse-field Ising chain, and a dis- 
cussion of its properties occupies the remainder of this 
subsection. A great deal is already known about the 
transverse-field Ising chain, and the following discussion 
is essentially a summary of known results as they relate 
to the KLM. Of particular importance, it is known that 
ff cr it undergoes a quantum phase transition from a FM 
phase, to a disordered PM phase. We will outline how 
the transition is formally determined below, but before 
proceeding it is perhaps useful to consider the physics of 
the FM-PM transition in the KLM. 

-ff cr it describes the double-exchange FM ordering be- 
ing gradually destroyed as the coupling J is lowered. 
The first term of ff C rit, which describes double-exchange, 
has been discussed extensively in subsection IIB above. 
The destruction of the double-exchange ordering may be 
understood physically as follows: As the coupling J is 
decreased the conduction electrons become less strongly 
bound to the localized spins, and tend to extend over 
spatial ranges beyond the effective range a for double- 
exchange ordering. Double-exchange becomes less ef- 
fective, and regions of ordered localized spins begin to 
interfere as the conduction electrons extend. The in- 
terference leads to spin-flip processes, and are embod- 
ied in -ff cr it in the transverse-field (the second term of 
Eq. (B2fl ) . The transverse-field in ff cr it includes two low- 
energy spin-flip processes by which the conduction elec- 
trons disorder the localized spins. One spin-flip process 
is backscattering, and is accompanied by a momentum 
transfer of 2kp from the conduction electrons to the local- 
ized spins. Since the chain of localized spins will tend to 
order so as to reflect this transfer, the transverse-field cor- 
responding to backscattering spin-flips is sinusoidal with 
modulation 2kp. The other low-energy spin-flip process 
in ff cr it is forward scattering. This involves zero momen- 
tum transfer to the localized spins, and the corresponding 
transverse- field is a constant (i.e. has modulation zero). 

It will become clear below that either forward or 



backscattering spin-flip processes separately are sufficient 
to destroy the FM order, and bring on a FM-PM phase 
transition in the KLM. However, for incommensurate 
conduction band filling, the backscattering spin-flip pro- 
cesses introduce a competing periodicity in the chain of 
localized spins. It turns out that this has non-trivial 
consequences for certain properties of the localized spins 
near the transition. In the following we first indicate 
how the FM-PM transition is determined for arbitrary 
transverse-fields. We then compare the properties of the 
localized spins which are disordered due to forward scat- 
tering (constant transverse-field), with the properties in 
which the spin disorder is due to backscattering (incom- 
mensurately modulated transverse-field) . We point out 
that the special properties resulting from incommensu- 
rate backscattering are at least qualitatively reproduced 
by treating the full transverse-field in ff cr it as a random 
variable with the appropriate (displaced cosine) distribu- 
tion. The subsection concludes with a summary of the 
properties of the random transverse-field Ising chain as 
they relate to the transition region in the KLM. 



Determination of the phase transition 

Using a Jordan- Wigner transformation to spinless 
fermions a^alj, and in the thermodynamic limit, ff cr it 
may be written 



ffcrit = { a ] A Jl a l + \ { a ] B 3l4 + h -' 
3,1 L 



(23) 



to an additive constant, where (Aji) and (Bji) are 
real symmetric and antisymmetric matrices, respectively, 
with non-zero entries 

Ajj = hj = A— {I + cos(2k F ja)} 
Ajj+i = = Bjj+i = -Bj +1 j = -J/4. 

The quadratic form Eq. ( p3] ) may be diagonalized for any 
transverse- field hj by using the method of Lieb et al. [ p6| . 
This gives 



ffcrit — 



E 



(24) 



to an additive constant, where rjt, rj k are creation and an- 
nihilation operators for free spinless fermions, and where 
the energies ujj, are eigenvalues of the symmetric matrix 
(A+B)(A — B). As the coupling J is decreased, ff cr it un- 
dergoes a quantum order-disorder transition from a FM 
phase, to a quantum disordered PM phase, signalled by 
the breakdown of long-range correlations between the lo- 
calized spins, and a continuously vanishing spontaneous 
magnetization. The critical line for the transition is de- 
termined by the critical coupling J c which solves [E7| 
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N 



J N ~2 N ~[[ ^ = 



(25) 



as N — > oo. The free energy of the localized spins be- 
comes non-analytic at points satisfying Eq. (p5|). 

The FM-PM transition at the coupling J c is generic 
to transverse-field Ising chains, and does not assume a 
particular form for the transverse-field hj |27j]. For ex- 
ample, if we consider only forward scattering spin-flip 
processes in the KLM, the transverse- field hj = AJa/a 
is a constant. Solving Eq. (|2^) with hj = AJa/a gives 
the quantum critical line for the FM-PM transition at 



87r 2 Asin(7m/2) 
a f^°dk cos(fca)A 2 (fc) 



(26) 



A detailed discussion of the properties of the Ising chain 
with a constant transverse-field, which describes the FM- 
PM transition in the KLM with backscattering neglected, 
is given by Pfeuty As a second example, if we con- 
sider only backscattering spin-flip processes in the KLM, 
the transverse-field hj = AJacos(2kpja)/a is sinusoidal. 
In real heavy-fermion materials, the number of available 
conduction electrons per localized spin will in general be 
irrational |p9| . In this case the transverse-field due to 
backscattering has an incommensurate modulation 2kp 
with respect to the underlying lattice of localized spins. 
Nonetheless a FM-PM transition still occurs. As shown 
in Ref. [30], the solution of Eq. ( p5| ) for incommensu- 
rately modulated transverse-fields yields a coupling J c as 
in Eq. (^6|) for the constant transverse-field. Thus the 
critical line for the FM-PM transition in the KLM with 
only backscattering spin-flip processes coincides with the 
critical line for the transition in the KLM with only for- 
ward scattering. 



Effects of the form of the transverse-field 

While the FM-PM transition itself is largely indepen- 
dent of the details of the transverse-field, there are signif- 
icant differences in the properties of the localized spins 
on either side of the transition depending on the par- 
ticular form of the transverse-field. The differences are 
most clearly apparent in the wavefunctions correspond- 
ing to the free fermions r/k of the diagonalized hamilto- 
nian Eq. (|^). The wavefunctions are always extended, 
or Bloch-like, for a constant transverse-field pq] . For 
an incommensurately modulated transverse-field, the be- 
havior of the wavefunctions is far more complex. The 
Ising chain with an incommensurate transverse-field has 
been studied extensively by Satija et al. ]|(] p2^ . The 
model is important as it has localized states in ID, and 
thus provides a link to random systems | p2[ . Numerical 
studies pfH show that the wavefunctions corresponding 
to the free fermions r)k are localized in the disordered 



PM phase, and undergo a spectral transition at the FM- 
PM phase boundary. In the FM phase, the wavefunc- 
tions are self-similar and the eigenvalue spectrum forms 
a Cantor set. Since the correlation functions for the lo- 
calized spins are determined by the wavefunctions, the 
KLM with backscattering possesses far different proper- 
ties to the KLM with only forward scattering spin-flip 
processes, even though both undergo a FM-PM transi- 
tion. Differences in the eigenvalue spectrum Uk of the 
diagonalized hamiltonian Eq. ( p4| ) lead to a similar con- 
clusion regarding thermodynamic properties; since the 
KLM with incommensurate backscattering has a fractal 
eigenvalue spectrum, its thermodynamics are far differ- 
ent to the KLM with forward scattering, which has the 
more standard (cosine-type) eigenvalue spectrum p8|| . 

The situation becomes yet more complex when we 
consider all possible low-energy spin-flip processes avail- 
able to the conduction electrons. -ff C rit includes for- 
ward scattering with zero momentum transfer, and is 
represented by a constant trans verse-field. i? C rit also in- 
cludes backscattering with an incommensurate momen- 
tum transfer 2k p, and is represented by a 2k p sinusoidal 
transverse-field. H cr i t does not include spin-flip interac- 
tions with momentum transfers at higher harmonics of 
2kp: at 4fcp, 6fcp, and so on. The higher harmonics will 
arise in a bosonization treatment which includes non- 
linear corrections to the conduction electron dispersion 
relation These corrections are very weak compared 
with the forward and backscattering spin-flip processes, 
and it is usual to neglect them. However, the addition of 
(even weak) higher harmonics in 2k f to the transverse- 
field hj will greatly alter the solution of Eq. (|2^) . Instead 
of one solution, there now occur an infinite number of so- 
lutions to Eq. (|25|), and these occupy a finite region of 
the parameter space jjo). The series of solutions is re- 
flected in numerical studies on the Ising chain with a 
transverse-field containing more than one incommensu- 
rate harmonic |30| , pl| ' The region of spectral transitions 
becomes broadened, and the wavefunctions correspond- 
ing to the free fermions r/k of Eq. ( pi]) are observed to 
undergo a cascade of transitions between extended, criti- 
cal, and localized behavior. The transitions in the wave- 
functions occupy a finite region of the parameter space, 
and coincide with the solutions of Eq. ( |25| ) at which the 
free energy becomes non-analytic. While the region of 
spectral transitions becomes broadened, this is not the 
case for the magnetic transition. The FM-PM transition, 
signalled by the vanishing of longe-range correlations be- 
tween the localized spins, is observed to remain sharp 

& 

The behavior observed by Satija et al. in the numer- 
ical studies discussed above is qualitatively identical to 
the behavior of the Ising chain with a random transverse- 
field. (Properties of the random transverse-field Ising 
chain are disussed extensively by Fisher p3fl , and are 
summarized below.) To see this identification, note that 
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the central feature of the random transverse-field Ising 
chain is that dilute regions of FM order may survive into 
the PM phase, and similarly that dilute regions of disor- 
der may continue into the FM phase. This feature is at 
the heart of Fisher's results Q, and is shared by H clit 
for incommensurate kp: As discussed above, a broad- 
ened region of spectral transitions about the true FM-PM 
transition occurs in the KLM with incommensurate con- 
duction band filling. Thus there are small regions in the 
PM phase in which the localized spins exhibit behavior 
normally associated with the FM phase, and vice versa. 
To further pursue the identification between a random 
transverse-field, and that present in H CI i t for the KLM, 
recall that the spectral transitions occur at points sat- 
isfying Eq. ( ^5| ) at which the free energy becomes non- 
analytic. There is an immediate identification between 
these non-analytic points, and the Griffiths singularities 
fjij present in random models, in which thermodynamic 
quantities such as the magnetization become singular in 
a range of parameter space about the non-random transi- 
tion. (See Ref. [33] for the Griffiths regions in the random 
transverse-field Ising chain.) 

The behavior of Zf C rit for incommensurate conduction 
band filling admits of a natural physical interpretation. 
The conduction band does not share the periodicity of 
the lattice of localized spins, and is unable to either 
totally order or totally disorder the lattice as the FM- 
PM transition is crossed. There remain dilute regions 
of double-exchange ordered localized spins into the PM 
phase as only a quasi-commensurate fraction of the con- 
duction electrons become weakly-bound, and become free 
to scatter along the chain, at the FM-PM transition. The 
remaining ordered regions are dilute enough that no long 
range correlations remain, but their existence dominates 
the low-energy properties of the localized spins near the 
transition. 

These considerations lead us to treat the transverse- 
field hj — AJa{l + cos(2kpja)}/a of -ff C rit as a random 
variable, so that hj is chosen from the displaced cosine 
distribution p(h)dh where 



p(h) 



1 



nAJa y'l - (ah/AJa — l) 2 



(27) 



As discussed above, this treatment of hj does not alter 
the basic FM-PM transition described by H cr a, and thus 
is not needed in order to plot the phase diagrams of the 
KLM in Figs. 3 and 6 below. However, it does account 
for the properties of the localized spins near the FM-PM 
transition as observed by Satija et al. in their numerical 
simulations. We conclude by noting that at low conduc- 
tion band filling the treatment of hj in Eq. ( |27j ) follows 
an analogous treatment in spin glass systems (cf. Ref. 
[3])- 



Properties of the localized spins near criticality 

Results on the random transverse-field Ising chain may 
be obtained from Fisher [ [33"| , who uses an approximate 
real-space renormalization-group (RG) analysis, which 
nonetheless yields asymptotically exact results at low- 
temperatures near criticality. Following Ref. [33], the 
critical coupling for the KLM is given by 



Jc 



47r 2 Asin(7rn/2) 
a J °° dk cos(fca)A 2 (fc) 



(28) 



The critical line thus retains the form of Eq. ( |26| ) for for- 
ward or backscattering separately, but is down by a factor 
of 2 as both spin-flip disorder processes are included. It 
is convenient to measure deviations from criticality by 

il 



5 = [var(log h)] 1 log ■ 



47r 2 Asin(7rn/2) 
J a J °° dk cos(fca)A 2 (fc) 



(29) 



where the measure of randomness is 



2n-l 



„ , \ 1 1.3. ■• -.(271-1) V- 1 i 2. 

var(logM = E { ~ 2 ,4....^ £ ™ \~ lo § 2 ' 



n=l 



(2n) 



8 = on the critical line, is positive in the disordered 
PM phase, and negative in the FM phase. 

The distinctive feature of Fisher's RG analysis is that 
it focuses on anomalous clusters of double-exchange or- 
dered localized spins which survive for small 5 into the 
PM phase, and similarly, rare disordered regions close to 
criticality in the FM phase. These are due to the incom- 
mensurability of the conduction band filling with respect 
to the lattice of localized spins, and the consequent in- 
ability of the conduction band, as a single many-body 
entity, to either totally order or totally disorder the lo- 
calized spins as the transition is crossed. It is the anoma- 
lous ordered (disordered) regions of localized spins in the 
PM (FM) phase which are responsible for the Griffiths 
singularities. Although these anomalous regions are very 
dilute, they dominate the low-energy properties of the 
spin chain. Thus, while typical correlations are much as 
in the constant transverse-field Ising chain, the measur- 
able mean correlations are dominated by the anomalous 
regions, and consequently greatly alter the low-energy 
behavior. 

An important prediction of our theory of the phase 
transition in the KLM is that the spontaneous magne- 
tization grows continuously from criticality into the FM 
phase: For the random transverse-field |53| 



Mo (J) ~ (-5) 







5<Q. 



(30) 



where (3 = (3 - y/E)/2 « 0.38 This disagrees with 
numerical diagonalization results on small systems ||, 
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which see a discontinuous jump in Mq at least at larger 
fillings, but note that regions of intermediate Mq have 
been observed in related studies and in small sys- 
tems at lower fillings 0. Indeed a discontinuous jump in 
Mq immediately above the transition seems difficult to 
understand in a thermodynamically large system, given 
that the ordering is due to double-exchange, and that the 
electron spin degrees of freedom are not frozen until deep 
into the FM phase |,|6|. 

Using Ref. [33], we summarize the properties of the ID 
KLM which are relevant to the transition region of small 
S. The mean spin-spin correlation function is defined by 



C{x) — (SjSj +x ) , 

where the average is over p(h), and where for convenience 
x denotes a continuous and positive variable. C(x) is 
dominated by atypically large correlations and for small 
| 5 | in the FM phase decays as 

C(x) ~ M 2 {5) +const.|<S| 2 ^/:r) 5 / 6 e~ 3 ^/«) 1/ V K/ « 

for x ^> £. The correlation length £ w 1/5 2 . (For typical 
pairs of spins the correlation length £ w S^ 1 ; the expo- 
nent is the same as in the Ising chain with a constant 
transverse- field.) At criticality, 8 = 0, the decay is power 
law: C(x) ~ x~P as x — ► oo. For small 5 in the PM 
phase, 



(31) 



where x ^> £ w l/<5 2 . Note that C (a;) decays more rapidly 
to Mq(S) in the FM phase, than it decays to zero in the 
PM phase. 

At low temperatures T close to the transition, C(x,T) 
decays exponentially at large distances with a correlation 
length £r. In the FM phase, £r diverges as a continuously 
variable power law of T: 



6 



o 2r r|i5 



/AS 2 T T \S\ 



oo 



in the FM phase, where Tt is a characteristic scale, given 
by Tt = log(max{ l /, hj}/T) at fixed J and n close to the 
transition. At criticality the correlation length is £t ~ 
AT^/tt 2 , while in the PM phase 



(6 2 + n/T 2 T y 



r T <5> i 



The correlation lengths for the long-range exponential 
decay of the correlations C(x,H) in small applied fields 
H along z have identical functional forms to those of £t 
above. We note only that in the FM phase, ~ H~ 2 \ s \ 
as H — > 0. This reflects the development of long-range 
order, and shows a power law dependence on H. (See 
Fisher ]3j| for more details). 

The magnetization in small positive applied fields H 
along the z direction is obtained |33| using an exact crit- 
ical scaling function. Close to the critical line in the FM 
phase this gives 



M(8,H) ~ M (tf)[l + 0(SH 2 ^ \ogH)} 

at zero temperature. At criticality, M(S, H) 
| \ogH\~P, and in the PM phase 



(32) 



Close to the transition in both phases the magnetization 
is highly singular. In the PM phase the magnetization 
has a power law singularity with a continuously variable 
exponent 26, and the linear susceptibility is infinite for a 
range of 5 into the PM phase. The susceptibility remains 
infinite (with a continuously variable exponent) close to 
the transition into the FM phase. The low temperature 
linear susceptibility x{T) takes the form 

in the FM phase, and T\{T) diverges as T -> 0. Note 
that the latter property was conjectured by Troyer and 
Wiirtz on the basis of their quantum Monte Carlo results 
@. At criticality, X (T) ~ T~ 1 \ logT| 2 ( 1_/3 \ and close to 
the transition in the PM phase, 



X (T) ~ S-^-^-^logT) 2 . 



(33) 



Tx(T) vanishes as T in the PM phase. The zero-field 
specific heat at low temperatures close to the transition 
is given by 

C V (T) - \S\ 3 T 2 ^[1 + 0{T) 2 ^} 

in either the FM or PM phases. At criticality, the specific 
heat C V (T) - |logT|- 3 . 



B. Weak-coupling 

Well below the transition, where Ja/vp is small and 
the FM interaction Eq. ( |l3| ) is negligible, the ordering 
of the localized spins is governed by the last two terms 
of the effective hamiltonian Eq. (^0|). To determine the 
dominant correlations in this strongly disordered phase, 
it suffices to take eigenvalues for ei(j) in the long-range 
object K(j) (cf. Eq. (21)). K(J) then fluctuates about 
zero incoherently, depending on the global Sj configura- 
tion. The effective hamiltonian corresponds to free lo- 
calized spins in x and z fields determined by conduction 
electron scattering. The free localized spin problem is 
straightforwardly diagonalized by standard methods jj7| , 
and yields a ground-state Sj configuration \ipo) given by 



exp 



i tan 1 



cos [if (j)] + cos(2fcpja) 
sin[if(j)] sin(2fcF ja) 



I) 



where | J.) is the state with S| 



-1/2 for all j. The 



dominant 2kp modulations in \ipo) are manifest, and are 
superimposed on an incoherent backgound: 



11 



(ii \S*S* +x \4> ) « sm[2k F ja] sm[2k F (ja + x)} 

to an incoherent normalization. This is observed at weak- 
coupling in numerical simulations ^|||[7j , and is called the 
RKKY regime. (Note, however, that the RKKY interac- 
tion strictly diverges in ID, and there is no lower bound 
on the ground-state energy for the RKKY hamiltonian, 
even for arbitrarily small J ]r2| . The divergence is typi- 
cal of perturbation expansions in ID, and does not occur 
in higher dimensions.) 



C. Phase diagram 

The behavior identified in the previous subsections is 
in complete qualitative agreement with the results of nu- 
merical simulations on larger systems ^fH?]] . To estab- 
lish quantitative agreement, i.e. to plot the critical line, 
we are presented with two obstacles. The critical line 
Eq. (H|) with both forward and backscattering spin-flip 
interactions included may be written 



J c a ■■ 



2ir 2 Av F 



: J °° dk cos(fca)A 2 (fc) 



(34) 



The first obstacle in using Eq. ( |34| ) is somewhat triv- 
ial, and relates to the global scaling of the critical line: 
The number A comes from the normalization of the Bose 
representations for spin-flip and backscattering electron 
interactions. It depends significantly on the cut-off func- 
tion A a (fc), and moreover relates to the normalization 
of Bose representations only in the limit of long wave- 
lengths (cf. Eq. (0)). The second obstacle is the depen- 
dence a = a(n, J) which measures the effective range 
of the double-exchange interaction. This is a non-trivial 
quantity in a thermodynamically large system. In our 
previous work j|] , we made the approximation of neglect- 
ing any functional dependence of a on n and J, and de- 
termined A by a fit to numerically determined points. 
The resulting phase diagram || correctly gives the gen- 
eral ground-state phase diagram of the KLM for J > 0, 
and indicates schematically the regions where Griffiths 
singularities occur for incommensurate filling, together 
with the crossover to the strongly-disordered (RKKY- 
like) regime at weak-coupling. Here we provide a more 
detailed analysis, and use numerically determined phase 
transition points to determine the functional dependence 
of a on the critical line. Note that contrary to all pre- 
vious results, the results of this subsection rely crucially 
on numerics. 

We estimate first the constant A. From bosonization, 
we know that a > 0(k F )~ 1 , and so a will diverge as the 
filling n — ► 0. From Eq. (34|) it follows that 



Jr. 



2^ An , as n — > . 



(35) 



Note the agreement with the exact solution of Sigrist et 
al. [pdl for the KLM with one conduction electron; the 



system is FM for all finite J. Recall also that the exact 
solution gives a/a = y/2/J for small J. It follows that a 
diverges at criticality in agreement with the result from 
bosonization as n — > 0. Using numerical results for J c 
at the smallest available filling (J c = 0.455 at n = 0.2 
from the infinite-size DMRG simulation of Caprara and 
Rosengren j?]]), we conclude from Eq. ( p5| ) that 27r 2 A « 
0.7. 
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FIG. 2. Plot of the dimensionless parameter Jcl/vf, which 
characterizes double-exchange FM, against coupling J > 
for numerically determined FM-PM transition points J c : the 
filled diamond is the quantum Monte Carlo result for systems 
up to 24 sites from Ref. [4] ; open circles and squares are exact 
numerical diagonalization results for the 8 and 9 site chain, 
respectively, from Ref. [5]; the filled square is the DMRG 
result for systems up to 75 sites from Ref. [6] ; the filled circles 
are infinite-size DMRG results from Ref. [7]. J c cl/vf ~ 0.7 
for vanishing J. The straight line of best fit is given, and 
shows good agreement with the spread of numerical results, 
together with the expected result as J — > 0. 

We now account for the functional dependence of a, 
and plot the critical line, a enters the critical line equa- 
tion in the denominator of the right hand side of Eq. 
(|34|). This factor may be determined independently of a 
choice for the cut-off function A a (k) by using numerically 
determined FM-PM transition points. In Fig. 2 we plot 
the dimensionless parameter Ja/vp against J for avail- 
able numerical data Ja/vp characterizes double- 
exchange in our theory, and gives the denominator of Eq. 
( p4] ) at criticality. The functional dependence is linear, 
and we give in Fig. 2 the straight line of best fit. This 
line gives 27r 2 A = 0.65 as n — > 0, in agreement with the 
estimate « 0.7 from n — 0.2 given in the previous para- 
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graph. It is reasonable to conclude that the deviations 
in the numerically determined points for Ja/vF from the 
straight line are reflections of the different critical values 
determined in different simulations. The line of Fig. 2, 
together with Eq. determines the critical line at 



Jr. = 



1.3 sin(7rn/2) 



J > 0. 



1 - 0.6sin(7m/2) 
The resulting phase diagram is given in Fig. 3. 



(36) 
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FIG. 4. The effective range a of the double-exchange 
interaction in units of the lattice spacing against filling 
n on the critical line. An exponential cut-off function 
A a (fc) = exp — (a|fc|/2) has been chosen. The vanishing of 
the range at half-filling in the J < KLM leads to a diver- 
gence in the critical line as half filling is approached. 



IV. THE FM-PM TRANSITION IN RELATED 
MODELS 



FIG. 3. Ground-state phase diagram of the ID KLM with 
J > 0. The critical line is from Eq. (^), and uses the line of 
Fig. 2. Numerically determined critical points are as in Fig. 2. 
At incommensurate fillings, there are Griffiths singularities in 
the free energy in a finite region of the parameter space about 
the critical line. At small Ja/vF in the paramagnetic phase, 
the system presents an RKKY-like behavior with dominant 
correlations in the localized spins at 2kF of the conduction 
band. 

The line of Fig. 2 determines the effective range a 
of the double-exchange interaction on the transition 
line. Choosing the exponential cut-off function A a (k) = 
e -a\k\/2 £ or s i m pii c ity ; we have 

or^coU)Asw = 1+1 ° /Q)2 - (37) 

For this choice of cut-off function, the line of Fig. 2 gives 
a/a = \/2.l/ J at the transition. This compares with the 
result y/2/ J obtained in the exact solution of the KLM 
with one conduction electron |ll|] just above the critical 
point at vanishing J. The filling dependence of a may 
be determined by using Eq. (H) to write a/ a = y/2.1/J 
in terms of n. The result is plotted in Fig. 4. 



In this section we consider some variants of the ID 
KLM of Eq. (|l|), and use our methods to describe the 
FM-PM transition at partial conduction band filling in 
these models. We consider firstly the effects of inter- 
actions between the electrons, and show that for repul- 
sive interactions the phase boundary is always pushed 
to lower J values, in agreement with available exact and 
numerical results p4|p^ ]. Second, we consider the KLM 
with a FM J < coupling. Our method extends to this 
model, and predicts the same class of phase transition as 
for J > 0. Using available numerical results ||, we follow 
the analysis of Sec. Ill to determine the critical line and 
give the resulting phase diagram. 

A. Interacting conduction band 

To determine the effects on the ordering of the localized 
spins due to interactions between the electrons, consider 
adding to the standard KLM hamiltonian of Eq. (|l|) the 
Hubbard interaction term 

v = uY^ n ii n n ■ 

3 



13 



In terms of density operators, the interaction may be 
written 

V = W + P-(k))(P + (-k) + P-(-k)) 

k 

-(a + (k)+a_(k))(a + (-k) + a_(-k))}, 

where the charge and spin density operators are de- 
fined in Eq. (||). It will be sufficient for our pur- 
poses to consider only forward scattering contributions 
to V. (For weak repulsive interactions, U small and 
positive, backscattering interactions renormalize to zero 
p8| .) Within the Bose description, this is equivalent to 
attaching the weight A Q (fc) to the density fluctuations in 
V. The interaction then reduces to standard Tomonaga- 
Luttinger-type, with forward scattering interactions de- 
scribed by bosonic density operators. The pure conduc- 
tion band part Hq + V of the interacting KLM may now 
be straightforwardly diagonalized via a Bogoliubov trans- 
formation exp(Ss) where Sb = J2v=p a an d 

S„ = — log — 

2L \vp 

x E \["+(k)u-(-k) u-{k)u + {-k)]K 2 a {k). 

k>0 

The charge and spin velocities are 



v p = Vpy/l + Ua/nvp , 
v a = VF\ji — Ua/irvp . 

By comparison with the exact Bethe ansatz solution, 
these velocities are correct to leading order in U. Cor- 
rections to the velocities at stronger couplings are given 
by Schulz fl39j . (Note that the spin velocity v a does not 
go complex as U increases, but smoothly goes to zero.) 
Under Sb the Bose fields transform as 



4>u(j) = \/ v f/v u 4> v (j), 
= \fv v /v F 6 v {j), 



(38) 



while 



v—p,a 



4tt 



to an additive constant. Under Sb, the bosonized KLM 
with interactions between the electrons now takes the 
same basic form as the original bosonized hamiltonian of 
Eq. (m. The only differences are that the first term in 
Eq. (pf) is replaced by ^„ H V1 and that the Bose fields in 
the remaining terms are replaced by their scaled forms 
as in Eq. ( |38| ) . Proceeding much as in the original U = 
problem, we choose the transformation exp(S) where 




(39) 



fv 



and obtain a transformed hamiltonian similar to Eq. (|1J 
The important difference is that the prefactor of the 
double-exchange FM term of Eq. (|l3|) is increased: 



J 2 a 2 



2„2 



J 2 a 



1 



Ait 2 vf 4,tt 2 vf 1 — Ua/irvF 



(40) 



Following exactly the analysis of Sec. II, we obtain an 
effective hamiltonian for the localized spins. This deter- 
mines a quantum order-disorder transition between FM 
and quantum disordered PM phases, with a critical cou- 
pling J C {U) which is down from the U — critical cou- 
pling by a factor 1 — Ua/nvF- For stronger interactions 
between the conduction electrons, the spin velocity of 
Schulz |39| should be used in the transformation Eq. (|3S|). 

The effect of a repulsive Hubbard interaction between 
the electrons is then as follows. Double-exchange is 
characterized by the enhanced dimensionless constant 
J(X\Jvf /v&- The FM phase becomes more robust, and 
the FM-PM phase boundary is pushed to lower values 
of J. We expect this on physical grounds: The repul- 
sion between the electrons tends to keep the regions of 
double-exchange ordered localized spins from interfering. 
Spin-flip disorder processes are thereby reduced. Our re- 
sult is consistent with the numerical work of Moukouri 
et al. jl5| on the KLM with t-J interacting electrons; re- 
duced critical couplings J c ~ 0.8,1,1.2 are determined 
for fillings n = 0.5,0.7,0.9, respectively. Moreover, since 
v a — > as U — > oo j39|, so that J c — > 0, our result coin- 
cides with the rigorous result of Yanagisawa and Harigaya 
fll4f for infinite repulsive electron interactions. 



B. The KLM with a FM coupling 

The KLM with a FM coupling J < is an effective 
model for CMR materials |2j . The localized spins in this 
model are three ti g Mn d-electrons, and have spin 3/2. 
The properties of the model of interest here are largely 
independent of the magnitude of the localized spins, and 
they may be approximated by spins 1/2 Noting that 
the double-exchange FM term Eq. ( |l3| ) is insensitive to 
the sign of J, it may be readily verified that the deriva- 
tion of the effective hamiltonian of Eq. ( p0| ) carries over 
to this case with minor modifications. This determines a 
quantum order-disorder transition from a FM to a disor- 
dered PM phase, with a critical line for \ J C \ given by Eq. 
©■ 

To plot the critical line, we follow the analysis of Sec. 
Ill and use available numerically determined transition 
points for the J < KLM. Yunoki et al. ^ determine 
the FM-PM transition for classical spins via Monte Carlo, 
and for quantum spins 3/2 via the DMRG. The resulting 
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transition lines, with coupling J correspondingly scaled, 
are very close, and their points may be used within our 
spin 1/2 approximation. In Fig. 5 we plot the dimen- 
sionless parameter Ja/vF against J for numerically de- 
termined points. The straight line of best fit gives very 
good agreement with the points, and as in Sec. Ill, de- 
termines the critical line 



Jc 



0.7sin(7rn/2) 
1 — sin(7m/2) 



(41) 



The resulting phase diagram is given in Fig. 6. The crit- 
ical line diverges close to half-filling, and differs from the 
J > KLM for which the line remains finite. We have not 
included the phase separated region identified by Yunoki 
et al. H in Fig. 6. Phase separation is observed in the 
classical spin simulation in the PM region from J c = 4. 
It is not observed in the quantum simulation until J = 6, 
and then occurs away from the FM-PM transition closer 
to half filling. Any phase separation involves strongly lo- 
calized electrons, and on-site localization is not described 
well by our bosonization of the conduction band (cf. Sec. 
II). 
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FIG. 5. Plot of the dimensionless parameter — Jcl/vf, 
which characterizes double-exchange ordering in the J < 
KLM, against coupling J for numerically determined FM-PM 
transition points: Open circles are results on classical local- 
ized spins using Monte Carlo on systems up to 40 sites from 
Ref. [8]; filled squares are DMRG results on a 16 site chain 
for quantum spins 3/2, and a correspondingly normalized cou- 
pling, from Ref. [8]. The straight line of best fit gives very 
good agreement with all points. 
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FIG. 6. Phase diagram for the ID KLM with a ferromag- 
netic coupling J < 0. The critical line is from Eq. did"]), and 
uses the line of Fig. 5. Numerically determined transition 
points are as in Fig. 5. Properties of the localized spins close 
to criticality and at weak-coupling are as for the J > KLM. 
The phase separated region identified in Ref. [8] for the clas- 
sical spins from J c — 4 into the paramagnetic phase is not 
shown. Phase separation is observed in Ref. [8] in the quan- 
tum simulation only at stronger couplings, and away from the 
FM-PM transition closer to half filling. 

The effective range a of the double-exchange inter- 
action on the transition line may be determined as in 
Sec. III. For an exponential cut-off function, we find 
a/ a = ^/0.7/ J. As a function of filling, this relation 
may be used together with the critical line Eq. (|4l]) to 
plot a I a against n as shown in Fig. 4. The vanishing of 
the effective range close to half filling is the reason the 
critical line diverges. 

The different filling dependence of a for J > and 
J < is shown in Fig. 4. Different effective ranges a for 
the double-exchange interaction for different signs of the 
coupling is due to the different infinite \J\ symmetries of 
the sites containing localized conduction electrons. For a 
FM coupling, this is a triplet with energy — | J|/4, while 
for Kondo couplings J > the on-site symmetry is sin- 
glet with the lower energy —3 J/4. The gain in energy for 
double-exchange per site and per conduction electron in 
either case is — | J|/4. For J < 0, the system thus gains 
just as much energy from double-exchange as it does by 
forming localized triplets, and the triplets have minimal 
effect for a non-vanishing hopping. The effective range 
a vanishes smoothly as the number of excess localized 
spins declines: a — > as n — > 1. For J > the situa- 
tion is far more complex. The on-site singlet energy is 
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lower than the gain for double-exchange, and there is a 
complicated co-existence between the two effects. (This 
is why the exact solution for the KLM with one conduc- 
tion electron is complex for J > 0, whereas it is trivial 
for J < Jul.) Fig. 4 indicates a saturation a « a as 
n — ► 1 for J > 0, in contrast to the J < behavior. 
We leave as an outstanding problem the reason for this: 
An accurate determination requires a detailed descrip- 
tion of localized Kondo singlet formation on a par with 
the double-exchange ordering and weak-coupling spin-flip 
scattering which have been the focus here. 

V. DISCUSSION AND CONCLUSIONS 

In this paper, we have described double-exchange FM 
ordering in the partially-filled ID KLM, and the destuc- 
tion of the FM phase by spin-flip disorder scattering. At 
weak-coupling deep in the disordered PM phase, the scat- 
tering determines RKKY-like correlations in the localized 
spins. Kondo singlet formation has been taken into ac- 
count indirectly, via an effective range for the double- 
exchange interaction. The effective range was identified 
as follows: A length a > 0(kp)~ 1 originates in bosoniza- 
tion as the minimum wavelength for density fluctuations 
which satisfy bosonic commutation relations. Bosoniza- 
tion describes fluctuations beyond a, and keeps the con- 
duction electrons finitely delocalized over lengths below 
a. The electrons preserve their spin over this range. We 
showed in Sec. II that this finite derealization may be 
identified with the length for coherent conduction elec- 
tron hopping, and measures the effective range of the 
double-exchange FM ordering induced on the localized 
spins by the electrons (cf. Eq. ( |l3| ) and Fig. 1). The rea- 
son this works is that double-exchange is conceptually 
a simple interaction. It reflects only the tendency for 
hopping electrons to preserve their spin, as they move 
to screen the more numerous localized spins. Double- 
exchange is characterized by the dimensionlcss factor 
Ja/vp, with Vp the conduction electron Fermi velocity 
and a the lattice spacing. 

We obtained a FM double-exchange interaction term 
Eq. (|l3| ) between the localized spins. The term was 
derived using a unitary transformation and is non- 
perturbative. This contrasts with other interactions de- 
rived for the localized spins in the KLM, such as the 
RKKY interaction, which are perturbative. The uni- 
tary transformation generates an effective hamiltonian 
Eq. ( p0| ) for the localized spins. The competing affects 
on the spin ordering are made manifest in the effective 
hamiltonian. The competing effects are double-exchange 
ordering at stronger coupling, and spin-flip disorder pro- 
cesses involving nearly free electrons at weak-coupling. 
The transition from a double-exchange ordered FM phase 
to a quantum disordered PM phase was then shown in 
Sec. Ill to be the quantum order-disorder transition of 



the transverse- field Ising chain (cf. Eq. (^6|)). This de- 
scribes double-exchange ordered regions of localized spins 
being destroyed as the electrons become weakly-bound, 
and become free to move and scatter along the chain. As 
the coupling J is lowered, the transition is signalled by 
a continuously vanishing spontaneous magnetization Eq. 
( |30| ) , and a breakdown in long-range correlations between 
the localized spins Eq. ([H]). The phase diagram is given 
in Fig. 3. Well below the critical line, no remnants of the 
ordering remain, and the effective hamiltonian describes 
dominant correlations in the localized spins at 2k p of the 
conduction band. 

Spin disorder occurs through forward and backscatter- 
ing spin-flip processes between the electrons and the lo- 
calized spins. We identified interesting properties result- 
ing from an incommensurate modulation of the backscat- 
tering momentum transfer with respect to the underlying 
lattice of localized spins: For incommensurate fillings, 
the conduction band has a competing periodicity with 
respect to the spin chain, and the electrons are unable to 
totally order, or totally disorder the spin chain at criti- 
cality. This leaves anomalous regions of double-exchange 
ordered localized moments close to criticality in the PM 
phase, as only a quasi-commensurate fraction of the elec- 
trons become weakly-bound at the transition. Similarly, 
there remain anomalous disordered regions close to criti- 
cality in the FM phase. The anomalous regions are very 
dilute, but they dominate the low-energy behavior of the 
localized spins. The magnetization Eq. ( |32| ) is highly 
singular for a finite range of couplings about the criti- 
cal line: The magnetization has a continuously variable 
power law exponent, and the susceptibility is infinite for 
a finite range of couplings even in the PM phase (cf. Eq. 

We considered in Sec. IV the effect of conduction elec- 
tron interactions on the FM-PM transition, and found 
that double-exchange is characterized by the dimension- 
less factor Ja^/vpjv^ for repulsive interactions, where 
v a is the conduction electron spin velocity. This factor 
is enhanced (cf. Eq. (fiol)) over the factor Ja/vp char- 
acterizing double-exchange with no interaction between 
the conduction electrons. This pushes the critical line to 
lower values of the coupling J, and for infinitely strong 
repulsive interactions FM occupies the entire phase di- 
agram for J ^ 0. The reason for this behavior is that 
for infinite repulsive interactions, the double-exchange 
ordered regions are prevented from interfering, and the 
spin-flip disorder processes arc ineffective. 

Since Kondo singlet formation is taken into account in 
our description only indirectly, our method extends also 
to the KLM with a FM J < coupling, and a FM-PM 
transition of the same class as J > was identified. The 
phase diagram for J < is given in Fig. 6. The difference 
between J > and J < is in the effective range a of 
the double-exchange interaction, and is due to the differ- 
ent energies for on-site triplets when J < 0, to on-site 
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Kondo singlets when J > (cf. Sec. IV). The effective 
range a, which enters bosonization as a finite but un- 
known length, was determined at the FM-PM transition 
by using our critical line Eq. (|34|), together with numer- 
ically determined transition points We found that 
a/a oc l/-\/]J], in agreement with a simple characteriza- 
tion of double-exchange at small conduction band fillings 
(cf. Sec. II) , and in agreement with an exact result jll| at 
vanishing filling. (Recall that the coupling J is measured 
in units of hopping t) . The proportionality constants are 
different for different signs of the coupling, and are fixed 
by a best fit to available numerically determined FM-PM 
transition points (Figs. 2 and 5). In Fig. 4 we plotted the 
corresponding filling dependence of a on the transition 
line: a — > as n — » 1 for J < 0, but remains finite for 
J > 0. This has a significant effect on the phase dia- 
grams. For J > (Fig. 3) the critical line remains finite 
as n — > 1, while for J < (Fig. 6) the critical line di- 
verges approaching half-filling. (Note that we do not con- 
sider the half-filled KLM. Indeed double-exchange FM is 
absent if the number of conduction electrons equals the 
number of localized spins.) 

The transition we identified is generic to partially-filled 
spin 1/2 KLMs, at least in ID. Our use of bosoniza- 
tion prevents us from anything more than speculation on 
the FM-PM transition in higher dimensional KLMs. We 
note only that (i) double-exchange is not restricted to 
ID f p2]|lg| ], and should be considered in any discussion 
of partially-filled KLMs in higher dimensions, (ii) Nu- 
merical work H on the KLM with a FM coupling does 
present a FM-PM transition in higher dimensions, which 
is very similar to the transition in the ID case. 

We conclude with a simple physical picture, suggested 
to us by our results, which underlies the generic ground- 
state transition. At small fillings in the FM phase, spin 
1/2 Kondo lattices form a gas of spin polarons, with each 
electron dressed by a cloud of ordered localized spins. 
The spatial extent of the polarization cloud is the effec- 
tive range a for the double-exchange interaction. For 
J > the localized spins tend to align opposite to the 
spin of the conduction electron. For J < they tend 
to align parallel to the electron spin. As the coupling is 
lowered, the polarization clouds gradually extend and be- 
gin to interfere. The interference causes spin-flip disorder 
processes, which eventually destroy the FM order: The 
spin-flip processes free the electrons from their clouds of 
polarized localized spins, and this signals the onset of 
the FM-PM phase transition. At couplings just below 
the transition in the PM phase, the electrons are nearly 
free, and move through the system. They scatter from 
the localized spins as they move, and the spin chain is 
disordered. At weak-coupling, the localized spins retain 
dominant correlations at 2k p of the conduction electrons, 
superimposed on an incoherent background. This reflects 
the momentum transferred from the conduction band to 
the spin chain in backscattering interactions, together 



with incoherent forward scattering. 
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